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.^ DISCLAIMER 

Although* It Is assumed that all the prbcautions have been 
taken to chccVbut this program thoroughly, no responsibility 
Is taken by the. originator of this programfor any erroneous 
results, miscboceptiorts, or mitrepresentatlons that may appear 
•In th\$ program. Furtnermore, no responsibility Is taken by 
Autonetfps Tn^ustriaK Products for jihe correct reoroductions of 
thispj^gramrNo Warranty, express or Jngalied, Is extended 
by me 06^ bi' ajipllc^^^ 
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U INTRODOCTION 

It is often desirous to have a medium of correspondence between two sets 
of data, empirical or otherwise, when there exists a degree of linearity 
between the functional representations of the data. This program cal- 
culates the linear function of one curve which best approximates a second 
curve by the least mean squares method. One curve X can be used to 
generate curves Y^ by fitting X over each Y^ in the best linear fashion* 
The program was principally designed to economize on time and non-linear 
equipment in analog conplter empirical curve simulation. 

2. METHOD . 

2,1 Let curve X and curve Y be functions, eii5>irical or otherwise, of the same 
variable 2. Choose N points from curve Y (the curve to be generated) at 
Z^ and N corresponding points fi-om curve X at the same Z^» Setting 
Ya = A + BX, the problem now resolves in finding the values of A and B 
which make the mean square value of Y - Y^^ a minimum. The mean square of 
Y ^ Ya = 





^ (Yi)^ -2aE Ii-2B £ YiXi 

^ T=l 1^1 

N N N 

+ A^ + 2AB r Xi + B^ 2 (Xi)^ = F (A,B) 
N. H . . 

F(A,B} has a minimum with respect to A at ^^ (mean square (Y-Yq^)) == 0. 

STa— 

* Hence, A = Y .- BX. • Substituting for A in equation I and minimizing we 
have . ^ ^ 

N N /\ f 

The percentage error at each point is computed as 
(P.E.)i = 100(Yi - Yai) 0< i^ N 
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3* RESTRICTIONS 

3.1 Components Required 

The photoreader and typewriter are required. 

3»2 Subroutines utilized, 

3.2.1 Floating Point Input AN-007.1 

3.2.2 Floating point to fixed point output AN-015.1 

3.2.3 Decimal output AN-016.0 

3.3 2^ N— 10 where N is the nymber of sets of points ffom both curves. It 
is obvious that N varies with the complexity of the curve being generated. 
Thus, in generating a strcdght line fiom a straight line, only two sets 
of points need to be chosen along the common variable Z. A more complex 
curve would require N to be larger for increased accuracy. (Notes 
Judicious selection of Z^'s will improve approximation; i.e., points 
should be chosen at slope intersections on curve 7). 

3.4 No point may be chosen which has a value of zero. 

3.5 If curve X is a straight line with slope zero, 3.3 becomes 2<^ N~ 10. 
Failure to meet this condition results in a possible overflow. 

4. . USAGE 

4.1 Conpiter set-up. 

4.1.1 Set senses-switch C down if heading is desired. 
4-«l*2 All other switches are in normal position. 

4.2 Sequence of manual operations. 

4.2.1 Load and verif^r program tape. The program will print out heading, if sense 
switch C is on, and halt waiting for input data. If the tape does not have 
L 0651 aO S at the end, press start 1 to achieve the same results. 

4.2.2 Set sense switch 3 down, and press start 1. The typewriter will carriage 
return and print "N:". Enter, on the typewriter, the nuniber of pairs of 
points chosen from curves X and Y 1^ typing an integer N, 2 ~ N '^ 10. 

4.2.3 Begin entering data. Enter N values, 2i£ N^ 10, from the Y curve first. 
followed by the corresponding N values from the X curve. The numbers may 
b@ integers, fractions, or mixed numbers. 
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4.3 Conyputation 

4*3*1 When all the data has been entered, depresss the start button to initiate 
computation. The values of A and B of the ecjaation Y = A + BX wi3.1 be 
typed out. The program vdll also tjype the difference between Y and the 
actual curve generated at each discrete point, along with the percent 
erroo: at that point. The average percent error is also determined. 

A*U Error codes and restart procedures. 

^•K*X If a correction is to be made in the data, put sense switch D down and 

press start 2. Enter first, on the typewriter, the number of correct 
pieces of data preceding the incorrect entry. Then re-enter the correct 
piece(s) of data. After corrections have been made, put sense switch D 
up, press error reset and start 2 to begin calculations. 

4..4«2 Computation Errors 

No errors are anticipated. 

Aw5 Extent of Storage Requirements. 

A. 5 Program Locations 0470g to 1200g 

1250g to la?^ 

OOOlg and 0002^ 

4. 5 #2 Subroutine Locations OOlOg to OAj67g 
4.5.3 Data Locations I2OO3 to I2473 

4.6 Output Format. 

4.6.1 The output format is the dame as that shown in the exaii^>le problem where: 

(1) Iffi (1:S i:S 10) is the difference at point i between the 
actual curve generated and the Y curve. 

(2) P.£. is the percent error at the corresponding point, 
(3). B( Average) id the average percent error. 
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5.0 EXAMPLE 
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CURVE GENERATION 
THE PURPOSE OF THrS PROGRAM IS TO GIVE THE BEST APPROX I MAT I ON OF ONE 
CURVE AS A LINEAR FUNCTION OF ANOTHER V/HEN BOTH CURVES ARE FUNCTIONS 
OF THE SArC VARIABLE. THE PERCENT ERROR OF APPROX 1 MAT I ON I S ALSO DETER- 
MIN^Do 

2 .0001 

2,3 .1 

2*6 
2.9 
2.9 

A: 1.999§0(?^;33^ 
B: 3.0008^9^17 

DYI5-. .000102 P.E.« .005125 

DY2: .000117 P-^E:, .005091 

DY3:. .000037 P*£> .001^08 

DY4: .000037 P.E. .00lii-08 

DY5:- .OOOOSif ?.^^^ .001514 

DYb:- > .0000l|4 P.E.- ^00\^\\ 

E(AVERAGE):- .0000 
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P.E, 

P.E. 

P.E. 

?.£.*^ 

P.E.- 
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6.0 
6.1 
6.1.1 



6.1.2 



6.2 
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CODING INFORMATION 

Constants and their locations (Octal)* 

Fixed Point 

Constant Location 

1060 at B18 In 0473 

Floating Point 

Constant Location 

100.0 0i^70 

2.0 0476 

0.0 126U 

1.0 ' 1262 

Timing Estimate. 

(16 +. ION) seconds, When N is the number of points chosen flrom the common 
variable of the X and Y curves. 



